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A  Lanchester  Model  of  Submarine  Attack  on  a 
Carrier   Battlegroup 

James  N.  Eagle 

Department  of  Operations  Research 

Naval  Postgraduate  School 

Monterey,  CA  93943 

Abstract 

A  Lanchester  model  is  developed  for  a  battlegroup  ASW 
engagement.   Two  variations  are  included.   In  the  first, 
long-range  missile  firing  submarines,  short-range  missile  or 
torpedo  firing  submarines,  and  submarines  firing  only 
torpedoes  distribute  their  attack  uniformly  over  battlegroup 
escort  ships  and  carriers.   In  the  second  variation,  the 
attack  is  concentrated  on  the  carriers. 


A    Lanchester    Model    of    Submarine    Attack    on    a 
Carrier     Battlegroup 

James  N.  Eagle 

Department  of  Operations  Research 

Naval  Postgraduate  School 

Monterey,  CA  93943 

I .  Introduction 

A  Lanchester  model  is  developed  for  a  battlegroup  ASW 
engagement.   The  model  was  developed  for  the  Center  for  Naval 
Warfare  Studies,  Naval  War  College,  Newport,  RI . 

II.  General  Description 

A  hypothetical  Blue  battlegroup  consists  of  two  types  of 
platforms:   carriers  (Z)  and  escort  ships  (A).   The 
battlegroup  is  being  attacked  by  three  types  of  Red 
submarines :   torpedo  shooters  (V) ,  torpedo  or  short-range 
missile  shooters  (C) ,  and  long-range  missile  shooters  (0) . 
When  platforms  engage,  either  or  both  can  be  destroyed. 

The  Red  submarines  regularly  return  to  a  Red  Base  for 
reloading  and  repair.   The  rate  of  return  is  proportional  to 
the  rate  at  which  the  submarines  deplete  their  weapons.   Red 
submarines  returning  to  the  Red  Base  pass  through  a  barrier 
patrolled  by  Blue  submarine  units  (S) .   Submarine  vs. 
submarine  engagements  occur  at  this  barrier.   Red  submarines 
must  again  pass  through  the  barrier  after  leaving  the  Red 
Base.   (See  Figure  1.) 

All  submarines  are  also  subject  to  area  attrition  by 
opposing  ASW  forces.   That  is,  in  each  time  period  a 
specified  fraction  of  the  remaining  submarine  units  are 
killed.   The  Blue  submarines  and  escorts  are  resupplied  at 
specified  rates.   All  other  units  are  not  resupplied. 


C^^)^, 


/ 


S  resupply   rate 


Red  Base 


©  © 


©  © 


Blue  SSN 
Barrier 


fl  resupply  rate 


Figure  1 .   Submarine  Attack  on  a  Battlegroup 


Table  1.  gives  the  names  of  the  kill  probabilities  used 
in  the  models  (kpl-kp21).   Table  2.  give  the  remaining 
variable  names . 

The  initial  values  for  the  state  variables  are  A(0), 
Z(0),  S(0),  Oi(0),  O2(0),  Vi(0),  V2<0),  Ci(0),  and  C2(0)  -  all 
of  which  must  be  specified.   Two  variables  are  used  for  each 
Red  submarine.   Oi(t),  V"i(t),  and  Ci(t)  are  the  numbers  of  Red 
submarine  units  attacking  the  battlegroup  at  time  t,  while 
02(t),  V2(t),  and  C2(t)  are  the  Red  submarines  at  the  Red 
Base  . 
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Table  1.   Probabilities 

of  Kill  Given  an  Engagement 

Rv:   Range  of  initial  detection  in  a  V/A  or  V/Z  engagement 
(nm)  . 

Rc :   Range  of  initial  detection  in  a  C/A  or  C/Z  engagement 
(nm)  . 

Vrei:   Average  relative  speed  for  random  search  between  Red 
submarines  and  Blue  surface  ships  (nm/hr) . 

area:   Area  of  battlegroup  (nm^) . 

a0,ac:   Engagement  rates  for  missile  attacks  by  0  and  C 
submarine  units  (engagements/day) . 

l0,lv,lc:   Fraction  of  0,  V,  or  C  submarine  units  which  deploy 
from  the  Red  Base  each  day  (1/day) . 

lg :   Fraction  of  S  submarine  units  leaving  the  Blue  SSN 
barrier  each  day  (1/day) . 

rsA, rsg:   Resupply  rates  for  A  and  S  units  (1/day) . 

dV/dct,ds:   Number  of  torpedo  attacks  per  patrol  that  the 
weapons  load  allows  for  V,  C,  or  S  units  (unitless)  . 

d0/dcm:   Number  of  missile  attacks  per  patrol  that  the  weapons 
load  allows  0  or  C  units  (unitless) . 

deltat :   Numerical  integration  time  step  (days) . 

n0,nv,nc,ns:   Fraction  of  alive  0,  V,  C,  or  S  units  killed  per 
day  by  opposing  area  ASW  forces  (1/day) . 

9or9c:      Tne  inverse  of  the  number  of  A  escort  ships  required 
to  reduce  the  probability  of  kill  by  an  0  or  C  missile  attack 

by  e_1  (1/ship) .   (A  large  g  means  the  AAW  missile  defense  is 
weak;  a  small  g  means  a  strong  AAW  missile  defense.) 

Robar^ Rvbar/  Rcbar  -      Detection  range  by  an  S  submarine  unit  on 
0,  V,  or  C  submarine  units  at  the  Blue  SSN  barrier  (nm) . 

Rsubavg:   The  average  of  Robar^Rvbar/  and  Rcbar- 

L:   Total  Blue  SSN  barrier  length  (nm) . 

vels,velt:   Searcher  and  transitor  speeds  for  the  barrier 
encounter  (nm/hr) . 

Table  2 .   Variable  Names 


III.    The  Lanchester  System  Models 

The  models  are  completely  described  by  expressing  the 
time  derivatives  of  each  of  the  nine  state  variables  as 
functions  of  the  state  variables.   Then  given  initial  values 
for  each  variable,  the  values  at  all  subsequent  times  can  be 
found  by  numerical  integration. 

Two  Lanchester  models  are  described.   The  first  models  a 
Red  submarine  attack  on  all  units  of  the  Blue  battlegroup. 
The  second  models  an  attack  only  on  the  carrier  units. 

A.   Model  1  -  Battlegroup  Attack 

Here  the  A  and  Z  units  of  the  battlegroup  are  attacked 
with  missiles  by  the  Oi  and  Ci  submarines,  and  attacked  with 
torpedoes  by  the  V"i  submarines.   The  defining  differential 
equations  follow: 


dA       2RV  Vrel  A  -gQA  _  , 

dt~  =  "   area   Vl  A  kp8  "  a7z~  (a°  kpl2  e     0l' 

-  X+z"  (a°  kp6  e~9°  Cl)  +  rSA 


Comments : 

1.  Escort  units  (A)  are  killed  by  random  search  interactions 
with  V"i  units  plus  missile  attacks  by  Oi  and  Ci  units. 

2.  The  terms  aQ  and  ac  are  engagement  (i.e.,  attack)  rates 
of  the  two  missile-shooting  platforms.   Constant  engagement 
rates  are  consistent  with  a  constant  cueing  rate  (from,  say, 
an  overhead  sensor) .   Attacks  are  uniformly  distributed  over 
the  A  and  Z  units.   The  Oi  and  Ci  platforms  are  either  unable 

or  unwilling  to  concentrate  fire  on  Z. 

3.  The  two  e  ^   terms  decrease  the  effectiveness  of  the 
missile  attacks  as  the  number  of  defending  A  units  increases. 


f  -  -  ^T^  Vi  Z   kpl8  -  Jj  (a0  kp20  ."**  ox) 

-gcA 


A+Z 


(ac  kpl6  e  ^c   Ci) 


Comments : 

1.  Z  units  are  killed  in  the  same  manner  as  A  units,  but  the 
kill  probabilities  are  different. 

2.  There  is  no  resupply  of  carriers. 


dS  I  .    2Rsuba 

dT=  "  min V'~l7s 


^     VK^T    J        (lCP21    +(1-kP21)/d^) 

( lo02    +    1VV2    +    1CC2    +   ao0i(l-kpll) /d0    +   acCi ( l-kp5) /dcm   + 

2RV   Vrf*i  ,  2RV   Vrei  ,  \ 

-  ViA(l-kp7)/dv      +     — Z^JL       ViZ(l-kpl7)  /dv    J 


area  area 
-    nsS    -    13S    +    rss 


Comments 

2R<31,Hairrr  /  /Ve  1  <A2 

is  the  probability  that  an 


"average"  Red  submarine  crossing  the  barrier  is  detected  by 
an  S  unit.   (Ref:   Search  and  Detection.  Alan  Washburn,  1981, 
p. 1-11.) 

2.  (kp21  +  (l-kp21) /ds)  is  the  probability  that  an  S  unit 

leaves  the  barrier,  given  that  the  S  unit  detects  a 
transitor.   The  S  unit  leaves  because  it  is  either  killed,  or 
it  is  not  killed  but  runs  out  of  torpedoes  during  the  attack . 

3.  lo02  +  1VV2  +  lcc2  are  the  rates  of  O2,  V2,  and  C2  units 
leaving  the  Red  Base. 

4.  ao0i (1-kpll) /dQ  is  the  rate  the  Oi  units  return  to  Red 
Base  after  running  out  of  missiles.   The  (1-kpll)  term  is 
necessary  since  a  returning  Oi  unit  must  both  run  out  of 
missiles  and  not  be  killed  in  the  counter-attack  following 
the  last  missile  salvo. 


5.  Similarly,  acCi (l-kp5) /dcm  is  the  rate  Ci  units  return  to 
Red  Base  after  running  out  of  missiles. 

2RV  Vrei 

6.  — T~ZZ —  ViA(l-kp7) /dv  is  the  returning  rate  of  Vi  units 

area 

to  Red  Base  after  their  last  interaction  with  an  A  unit . 

2RV  Vrei 

7.   — —  ViZ  (l-kpl7)  /dv  is  the  returning  rate  of  Vi  units 

area 

after  their  last  interaction  with  a  Z  unit. 

8.  nsS  is  the  rate  of  S  kills  by  Red  area  ASW  forces. 

9.  1SS  is  the  rate  of  S  leaving  the  barrier  due  to  causes 
other  than  being  killed. 

10.  rs$  is  the  resupply  rate  of  S  units. 


-£-  =  loqo02  "  ao0i  (kpll  +  (l-kpll)/d0)  -  no0i 
where  q0  =  1  -  kP22  minjl  .^ttE ^  1+(^f   } 


Comments : 

1.  qG  is  the  probability  of  an  O  unit  surviving  a  barrier 
transit  against  S  Blue  SSNs .   Note  that  survival  could  result 
from  no  detection  by  Blue  SSNs,  or,  if  detection,  an 
unsuccessful  attack. 

2 .  loqo°2  is  the  rate  at  which  Oi  units  arrive  at  the 
battlegroup  from  Red  Base. 

3.  ao0i  (kpll  +  (l-kpll)/d0)  is  the  rate  at  which  Oi  units 
leave  the  battlegroup  area  due  to  either  being  killed  after  a 
missile  launch  or  expending  all  their  missiles.   (Remember: 
aQ  is  the  Oi  engagement  rate.) 

4.  no0i  is  the  rate  Oi  units  are  killed  by  Blue  area  ASW. 


dY- 

dt     -    lvqvV2 
2RV   Vrei 
area 

where    qv   =    1 

V] 

2RV   Vrei 

area 

Z     (kpl7+ 
kp22   min 

ViA     (kp7+(l- 
(l-kpl7)/dv) 

•kp7)/dv 
-   nvVi 

) 

/•/=:,  i 

>eij 

Comments : 

1.  qv  is  the  probability  of  a  V  unit  surviving  a  barrier 
transit  against  S  Blue  SSNs . 

2.  — - — —  V]A  (kp7+ (l-kp7) /dv)  is  the  rate  Vi  units  leave 

the  battlegroup  area  by  being  killed  in  engagements  with  A 
units  or  by  running  out  of  torpedoes  in  such  an  engagement. 

2RV  Vr-ei      /  x 

3.  Similarly,  -~z —  ViZ  (  kpl7+ ( l-kpl7  ) /dv)  is  the  rare  V: 

units  leave  the  battlegroup  area  due  to  engagements  with  Z 
units . 


4 .   nvVi  is  the  rate  of  area  ASW  kills  of  Vi 


dC 

"££"  =  icqcC2 

where    qc   =    1 

acCi  {  kp5    + 
-    kp22    min")  1 

(l-kp5)  /d0)     -   ncCi 

.^rr  V'-ft; 

And  at  Red  Base,  we  have, 


02 

—  =   -lo02    +    ao0iqo(l-kpll) /d0 

d"/2 

2RV   Vrei                                               2RV   V^i 

'    -1VV2    +                        ViAqv(l-kp7)/dv  +              ~         ViZqv (l-kpl7) /dv 
area                *             r                       area           *                  * 

C2 

—  =    -1CC2    +   acCiqcd-kpl7) /dcm 

Comment : 

1.  c[o/  ^v/  and  3c  are  defined  as  before. 

2.  The  rate  of  change  of  each  of  the  O2,  Vo,  and  C2  groups  is 
the  sum  of  the  outgoing  rate  (e.g.,  -10C>2)  and  the  incoming 
rate  due  to  0]_,  V]_,  or  Ci  running  out  of  weapons  and  then 
surviving  the  Blue  barrier. 

B .   Model  2 :   Carrier  Attack 

Now  we  postulate  that  the  attack  is  concentrated  on  the 
remaining  Z  units.   We  require  the  Ci  platforms  shoot  only 
torpedoes.   When  a  torpedo  shooter  encounters  an  A  unit,  it 
attempts  to  evade  and  break  contact.   Finally,  we  assume  that 
the  Oi  units  can  direct  all  of  their  fire  at  the  Z  units.   If 
this  is  not  desired,  the  model  is  easily  modified  to  allow 
any  probabilistic  split  between  A  and  Z. 


dA 
dt 


2RV  Vrei 
area 


V]A  kplO 


2RC  Vrel 
area 


CiA  kp4  +  rsA 


Comment : 

1.   The  A  units  now  receive  no  missile  fire;  and  if  kplO  and 
kp4  are  small  (as  would  be  expected) ,  then  A  will  suffer  very 
little  attrition. 


dZ 
dt 


2Ry  Vrel 
area 


ViZ  kp!8  - 


2RC  Vrei  -a^A 

— ^^   CiZ  kpl4  -  aQ  kp20  e  9o   0X 

area        r 


Comments : 

1.   The  V]_  and  Ci  units  are  conducting  random  search  for  Z 

units . 


All  of  Oi's  fire  is  directed  at  the  Z  units 


dS 

'it: 


J.     2R3ubava       ^  /"TTEIIIf"  1      ( y,    o 


■^.^Sr  V1+l™  i  (kp21  +(1-kp21)/ds) 

( lo02  +  1VV2  +  1CC2  +  ao0i (1-kpll) /d0  + 

2RV  VrPi  2RC  Vrei  \ 

E£i  ViZ(l-kpl7)/dv   + ^^  CiZ(l-kpl3) /dct  ) 

area  *  area  r      LL  ' 

-  nsS  -  1SS  +  rss 


Comments : 

1.  As  before,  S  units  are  killed  in  submarine  vs.  submarine 
encounters  with  transitors  going  both  ways  through  the 
barrier . 

2.  Transitors  going  to  the  battlegroup  do  so  at  rates 
proportional  to  their  numbers  at  Red  Base. 

3.  Transitors  going  to  Red  Base  do  so  at  rates  proportional 
to  the  rate  they  run  out  of  weapons. 


dOi  /  x 

-  =   loqo02    -   ao0i  (kpll    +    (1-kpll) /dQ)     -   no0i 


dt 


Comment : 

1.   Oi  units  are  now  aiming  at  Z  units  exclusively.  But  the 
rate  at  which  Oi  units  fire,  are  killed,  and  run  out  of 
weapons  remains  as  before.   So  there  is  no  change  in  the 
expression  for  dOi/dt . 


dVl         i        w  2RV  Vrei 

"dT  =    lvqvV2    "        area       VlA   kp9    " 

Rareael   VlZ     ( kpl7+ ( 1_kP17 } /dv)  "    nvvl 


Comment : 

1.   No  torpedoes  are  fired  by  Vi  units  in  V"i/A  interactions 


dCi  =  ,    n  2RC  Vrei 

dt 
2RC  Vrei 


^^2  "   area   VlA  kp3  " 


CiZ  (kpl3+(l-kpl3)/dct)  -  ncCi 

cll.  fc:  a. 
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Comment : 

1.   Both  Ci  and  Vi  units  are  exclusively  torpedo  shooters 

And  at  the  Red  Base, 


02_ 

dt 

=    - 

lo02    +    ao0iqo(l-kpll) /d0 

dV2 

dt 

=    - 

ivv2 

2RV   Vrei 
area 

•  ViZqv(l-kpl7) /dv 

£2. 

dt 

=    -1 

•CC2 

2RC   Vrei 
area 

CiZqc(l-kpl3)/dct 

Comment : 

1.   As  before,  the  total  rate  of  change  of  submarine  units  at 

Red  Base  is  the  arrival  rate  minus  the  departure  rate. 


11 


DISTRIBUTION  LIST 

NO.  OF  COPIES 

Library,  Code  0142  2 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

Director  of  Research  Administration  1 

Code  012 

Naval  Postgraduate  School 

Monterey,  CA  92943-5000 

Library,  Code  55  1 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

Defense  Technical  Information  Center  2 

Cameron  Station 
Alexandria,  VA  22314 

Commander  Submarine  Development  Squadron  Twelve  1 

US  Naval  Submarine  Base  New  London 
Groton,  CT  06340 

Dr.  Al  Roth  1 

Department  of  Economics 
University  of  Pittsburgh 
Pittsburgh,  PA  15260 

Center  for  Naval  Analyses  1 

4401  Ford  Ave. 
P.O.  Box  16268 
Alexandria,  VA  22302-0261 

Center  for  Naval  Analyses  1 

Attn:  John  Hall 
4401  Ford  Ave. 
P.O.  Box  16268 
Alexandria,  VA  22302-0261 

H.  A.  Miller  (Code  33T)  1 

Naval  War  College 
War  Gaming  Department 
Newport,  RI  02841-5010 

Dr.  Don  Daniel  1 

Center  for  Naval  Warfare  Studies 
Naval  War  College 
Newport,  RI  02841-5010 


Distribution  List   Continued  no.   OF  COPIES 

John   Hanley  \ 

Center  for  Naval   Warfare  Studies 
Naval   War  College 
Newport,  RI   02841-5010 

Dr.  Tony  Ciervo  1 

Pacific  -  Sierra  Research  Corp. 
12340  Santa  Monica  Blvd. 
Los  Angeles,  CA  90025 

Operations  Research  Ctr,  Rm  E40-164  1 

Massachusetts  Instructor  of  Technology 
Attn:  R.  C.  Larson  &  J.  F.  Shapiro 
Cambridge,  MA  02139 

Arthur  P.  Hurter,  Jr.  1 

Professor  and  Chairman 

Dept  of  Industrial  Engineering 

and  Management  Sciences 
Evanston,  IL  60201-9990 


3  2768  00333297  4 


